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ABSTRACT 

Gravitational waves from the coalescence of binary black holes carry away linear momentum, causing 
center of mass recoil. This "radiation rocket" effect has important implications for systems with 
escape speeds of order the recoil velocity. We revisit this problem using black hole perturbation 
theory, treating the binary as a test mass spiraling into a spinning hole. For extreme mass ratios 
{q = mxjmi <C 1) we compute the recoil for the slow inspiral epoch of binary coalescence very 
accurately; these results can be extrapolated to g ~ 0.4 with modest accuracy. Although the recoil 
from the final plunge contributes significantly to the final recoil, we are only able to make crude 
estimates of its magnitude. We find that the recoil can easily reach ^ 100 — 200km/s, but most likely 
does not exceed '-^ 500km/s. Though much lower than previous estimates, this recoil is large enough 
to have important astrophysical consequences. These include the ejection of black holes from globular 
clusters, dwarf galaxies, and high-redshift dark matter halos. 

Subject headings: black hole physics — gravitation — gravitational waves — galaxies: nuclei 



1. INTRODUCTION AND BACKGROUND 

Along with energy and angular momentum, gravita- 
tional waves (G Ws) carry linear momentum away from a 
radiating source llBonnor fc R,otenbergll96H iPeresllTOB^ 
iBekenstdnl Il973j) . Global conservation of momentum 
requires that the center of mass (COM) of the system 
recoil. This recoil is independent of the system's total 
mass. 

iFitchetd l)1983(l first computed GW recoil for binaries. 
He treated the members as non-spinning point masses 
(mi,TO2), the gravitational force as Newtonian, and in- 
cluded only the lowest GW multipoles needed for mo- 
mentum ejection. For circular orbits Fitchett's recoil is 
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where rterm is the orbital separation where GW emission 
terminates, q = mi/m2 < 1 is the mass ratio, and M = 
mi -|- 7712 is the total mass. The function f{q) = q^{l — 
q)/{l + q)^ has a maximum /max at g ~ 0.38, is zero for 
q = 1, and has the limit f{q) « q^ for q 1. 

Equation |^ tells us that in the coalescence of binary 
black holes (BHs) — where rtcrm can approach GM/c^ — 
the kick might reach thousands of km/s. This is far 
greater than the escape velocity of many globular clus- 
ters (typically ~ 30 km/s), and may even exceed galactic 
escape velocities 1000 km/s). Rec oil could thus have 
important astrophysical implications ijRedmount Sz Ree^ 
[some of which are explored in a companion paper 

Icrritt ct al. 2004'; Paper II)]. This has motivated us 
to revisit this problem. 

Equation indicates that the recoil is strongest at 
small separations, when the relativistic effects neglected 
by Fitchett are most important. This issue has been 
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address ed in restricted circumstances using perturbation 
theory (Nakamura & Haugan"'1983'; 'Fitchett & DetweileJ 
198 4: Nakamura, Oohara, & Koiima 1987), post- 
Newtonian expansions (jWiseman 1992; Kidder 1995), 
and numerical relativitv llAndrade fc PricCi J^^if 
Anninos fc Brandti Il998t iBrandt &: AnninosI Il999t 
Lousto fc Pric el l2004l) . Unlike previous analyses, our 
treatment applies to the strong-gravity, fast-motion 
regime around spinning holes undergoing binary coa- 
lescence. Using BH perturbation theory we model the 
dynamics of the binary, the generation of GWs, and the 
backreaction of those waves on the system up to the 
inner-most stable circular orbit (ISCO). Our results are 
accurate only for extreme mass ratio inspirals {q 1), 
but we can extrapolate to q ^ 0.4 with modest error. 
We model the GW emission from the final plunge more 
crudely. 

2. OVERVIEW OF GRAVITATIONAL RADIATION RECOIL 
The rate at which momentum is radiated is given by 

where /i+.x are the "plus" and "cross" GW polarizations, 
77'^' is a unit radial vector from the sour ce, and r is the 
distance to the observer l)Thornelll980 h. [We have set 
G = c = 1; an overdot refers to a derivative with respect 
to coordinate time t; angle brackets denote averaging 
over several wavelengths.] The binary's COM recoil is 



-dP^^/dt. 



"-^ GW 

dt 



Decomposing hj^^x into multipoles in the wave zone 
(jrhornc lESfl.), Eq. (j^J can be expanded (to lowest order) 
as 

2 / d'^T^^ d^T^ \ 16 / d^JP^ d^S'J^ 
63\~dt^ dt3~/^45y ~dt^ dfi' 

■(3) 

where X*^ , 5*^ , and X*-''^ are the symmetric, trace-free 
mass quadrupole, current quadrupole, and mass octupole 
moments. Recoil thus arises from "beating" between dif- 
ferent multipoles. Applying Eq. |j2Jl to a Newtonian bi- 
nary and integrating yields Eq. (Q. 
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Fig. 1. — Recoil from prograde coalescence with a/Af = 0.8, 
r] = 0.1 {q = 0.127). Solid (blue) lines represent quantities during 
the inspiral, as calculated using our Teukolsky equation solver. 
Dashed (red) lines are calculations during the plunge (using the 
"upper-limit" prescription discussed in section U). The plunge is 
truncated shortly before the particle enters the event horizon. The 
different panels are: (a) Orbit of the mass /i about the central 
spinning hole. The dashed circle is the location of the ISCO. (b) 
Recoil velocity of the center of mass. The spiral ends when GW 
emission is cut off. (c) Motion of the binary's center of mass, (d) 
Total center of mass recoil velocity, {v'^ + Vy)^^^ . 



iWisemanl I^1992^ provides an intuitive description of 
the recoil: When two non-spinning bodies are in circular 
orbit, the lighter mass moves faster and is more effective 
at "forward beaming" its radiation. Net momentum is 
ejected in the direction of the lighter mass's velocity, with 
opposing COM recoil. When mi = m2, the beaming is 
symmetric and the recoil vanishes. The instantaneous re- 
coil continually changes direction over a circular orbit, so 
the COM traces a circle. Neglecting radiation reaction, 
this circle closes, and the recoil averages to zero over each 
orbit. With radiative losses, the orbit does not close, and 
the recoil accumulates. This accumulation proceeds un- 
til the holes plunge and merge, shutting off the radiated 
momentum flux and yielding a net, non-zero kick velocity 
(cf. Fig.d. 

Spin complicates this picture by breaking the binary's 
symmetry. Consider an equal-mass binary, with one 
member spinning parallel to the orbital angular momen- 
tum. Due to spin-induced frame dragging, the non- 
spinning body's speed — and hence radiation beaming — is 
enhanced. Kidder ( 1995) has treated this spin-orbit in- 
teraction in post-Newtonian theory. Specializing his Eq. 
(3.31) to a circular, non-precessing orbit, the total kick 
for two bodies with spins Si^2 = 01,2^1 2^ parallel (or 
antiparallel) to the orbital angular momentum is 
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where the spin-orbit scaling function /so('?: Qii 'I2) 
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q {0.2 — qai) / {1 + q) . The "correction" causes significant 
recoil even when q — 1 (and hence Vf = 0). The spin- 
orbit term is largest when q ~ 1 and the spins are max- 
imal and antiparallel (ai = —02 = ±1; /so, max = 1/16). 
The recoil vanishes for q — I and spins equal and parallel 
(oi = 02) — a symmetric binary. 

Since we work in the g <C 1 limit, we ignore the smaller 
body's spin, which incurs an error ~ q^ai in the orbital 
dynamics (Kidder 1995). Our extreme mass ratio anal- 
ysis treats the binary in an effective-one-body sense: a 
non-spinning point particle with mass /i — 77111112/ M or- 
bits a Kerr hole with mass M = rrii + m2 and spin S = 
alVPz. There is an ambiguity, however, in how one trans- 
lates the physical spin param eter 0,2 of the h ole to the 
"effective" spin parameter a. iDamourl l)2001tl provides 
a relation between these parameters, valid in the post- 
Newtonian limit for d < 0.3: a = 02(1 4- 3q/4:)/{l + q^ . 
Because of this ambiguity, we present our results in terms 
of the effective-spin-parameter a. Even if the larger hole's 
spin is nearly maximal (5,2 — ±1), finite mass ratios 
q > 0.1 restrict our results to spins with \a\ < 0.8 — 0.9. 

When applied to a perturbation calculation of the 
head-on collision of two BHs, an effective-one-body scal- 
ing of the GW energy flux {Eqw oc q"^) in which q 
rj — /i/M — q/{\ + qY has been shown to agree with 
results from full numerical relativity ( Sm arr ..1978 ^. We 
use a similar "scaling up" procedure for the momentum 
flux: In perturbation theory P^^fj oc q^. We then sub- 
stitute q^ f{q) l^tchett _&Detwei^[lll3)- terms 
of rj, the scaling function is given by f{q) /{rj) — 
77^^/1 — 477, and is maximized at 77 = 1/5.] Using f{q) [or 
fifj)] to scale the momentum flux assumes both bodies 
arc non- spinning and that the orbit is quasi-circular. For 
simplicity, approximate spin corrections to /(g) based on 
Eq. Q are ignored (incurring errors < 30% if i? < 0.4) 
(cf. Paper H). 

3. INSPIRAL RECOIL FROM PERTURBATION THEORY 

Our model binary consists of a mass ^ in circular, equa- 
torial orbit about a BH with mass M and effective spin 
a — dM. (GWs rapidly reduce eccentricity, so circularity 
is a good assumption for many astrophysical binaries.) 
When n <^ M, binary e volution is well d escribed using 
BH perturbation theory ljTeukolskvlll973() . We treat the 
binary's spacetime as a Kerr BH plus corrections from 
solving the perturbed Einstein equations — the Teukol- 
sky equation. Specifically, we solve a linear wave equa- 
tion for the complex scalar function ^1/4, which describes 
radiative perturbations to the hole's curvature. Far from 
the binary, ^4 = {h+—ihx)/2; it therefore encodes infor- 
mation about the GW fields in the distant wave zone, as 
well as the energy, momentum, and angular momentum 
carried by those fields. 

Far from the binary 5*4 has the expansion: 



*4 = -^^i™5;„(0;ac^„)e™^ 



(5) 



Im 



In terms of Boyer-Lindquist coordinates (t, r, 9,(j)), = 
t — r in retarded time, uj^ = 7nrioib is a harmonic of the 
orbital frequency, Sim{0;aujm) is a spheroidal harmonic, 
and Z;„i is a complex number foun d by solving a partic- 
ular ordinary differential equation l)Hughesll2000() . 
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The linear momentum flux can be extracted by com- 
bining Eqs. (01 and (O- The resulting expression is sim- 
plest in the "corotating" frame, (ff™°'^ = (pit) — ftorht- 

■ 1 ^;m-Z^/'(m+l) 



ll'm 



'0 



(6) 

Here, T^gw — e *"^^*^[^gw + *^gw]' ^"^^ overbar de- 
notes complex conjugation. Similar expressions give the 
energy and angular momentum fluxes. The recoil veloc- 
ity is found by integrating Eq. © , starting at initial time 
To when the binary is at large separation [and the recoil 
is well described by Eq. iQJ] 
GW emission terminates: 



and ending at time T when 



1 

M 



Vgw dt 



(7) 



Our procedure starts with a point-source on a circu- 
lar geodesic orbit with specified energy E and angular 
momentum L^. Solving the Teukolsky equation gives us 
the energy, momentum, and angular momentum fluxes 
of GWs to infinity and down the event horizon. (The 
linear momentum flux down the horizon does not affect 
the recoil.) In the adiabatic limit (in which GW back- 
reaction changes the orbit very slowly, r/r <C 27r/r2orb), 
the energy and angular momentum fluxes (-Eqw- -^z,gw) 
are used to evolve to a new geodesic with E — EQ^^^At 
and Lz — Lz.Gw^t- Repeating this procedure for a se- 
quence of geodesies generates a slow inspiral trajectory. 
The momentum flux along this trajectory and associated 
recoil velocity are then calculated via Eqs. © and (0. 

This prescription can be used to calculate the recoil 
velocity only up to the ISCO. There the slow, adiabatic 
inspiral of the particle transitions to a rapid "plunge" 
that terminates when the particle crosses the event hori- 
zon (cf. Fig.^). Our Fourier decomposition of ^'4 is no 
longer valid as there are no well-defined harmonics lo„i 
for plunging trajectories. 

Figure shows the perturbation theory calculation 
of the ISCO recoil for a binary with reduced mass ratio 
77 = 0.1 {q = 0.127). The solid curve in Figure la can be 
fit by 



Msco = 422 km/s 



fma: 



2.63+0.06nsco/M 



(8) 



where risco is the spin-de pendent ISCO radius [defined 
for q = by Eq. 2.21 of iBardeen. Press, fc Teiikolskvl 
l)1972() ]. and we have included the appropriate scaling 
function (valid for q < 0.4). Although our adiabatic 
assumption is violated for rj = 0.1 (especially for large, 
prograde spins) our results are still valid since Visco / fiq) 
is only weakly dependent on q (and is independent of q 
in the g — *■ limit). 

For retrograde orbits around rapidly spinning holes, 
the ISCO is at large radius (9Af for a — —1) and Fitch- 
ett's Newtonian formula [Eq. ^] agrees well with our 
result. For prograde inspiral into rapidly spinning holes, 
the ISCO is deep in the strong field, where relativistic ef- 
fects become important and suppress the recoil relative 
to Fitchett's result. 

4. RECOIL ESTIMATES FROM THE FINAL PLUNGE 
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Fig. 2. — Recoil velocity versus effective spin a/M for ?j = 0.1 
((J = 0.127). (a) Recoil velocity up to the ISCO. The solid (blue) 
curve is our Teukolsky equation result. The dashed (red) curve 
shows the Newtonian recoil prediction [Eq. 0], which is sub- 
stantially higher for large, prograde spins (smaller ISCO radius), 
(b) Upper and lower limits for the total recoil. The shaded re- 
gion represents our uncertainty in the final kick velocity. The de- 
tailed shape of the upper-limit curve depends on the nature of our 
truncated-power-law ansatz. 



During the plunge, the small body's motion is dom- 
inated by the Kerr ef fective potential ra ther than 
radiation-reaction forces ijOri fc Thornel2000(l . It is easy 
to match a plunging geodesic with constant E and 
onto an inspiral trajectory near the ISCO. With a code 
that does not Fourier expand ^'4 (Khanna 2003; Mart^ 
l200l . one could properly compute the GW emission and 
associated recoil along such a plunging trajectory (when 
9« 1). 

Since we do not have such a code at hand, we must 
estimate the wave emission more crudely. Our results 
from the inspiral show that, for a given spin, T^gw is 
well described by a power law in radius, Pgw oc 
from large r up to the ISCO. As an approximate "upper 
limit" of the recoil, we make the ansatz that this power 
law can be continued past the ISCO. This must break 
down at some point: the power-law reflects the circu- 
larity of the inspiral orbit and should be suppressed by 
the increasingly radial motion during the plunge. To pre- 
vent the momentum flux from diverging, we truncate the 
power law at 3Af, replacing it with the condition that 
(dtj dryPiGSN ~ constant, where r is proper time along 
the plunge geodesic. This allows the momentum flux to 
"redshift away" as the particle approaches the horizon. 
Using the recoil velocity at the ISCO as initial condi- 
tions [Sec. 01 and a plunge trajectory with coordinates 
[r(t),0(t)], we use Eq. and our truncated-power-law 
ansatz to compute the accumulated recoil until a cutoff 
time T when the horizons of the holes come into contact 
(in a quasi-Newtonian interpretation of the coordinates). 
The upper-curve of Figure|2b shows the result of this cal- 
culation (for r\ = 0.1). 
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We also perform a separate "lower-liniit" calculation. 
A plunge trajectory is computed as before, but in place 
of the power-law ansatz for VqWi we integrate the trun- 
cated, multipole expansion of Eq. Q instead. In this cal- 
culation the momentum flux initially grows like a power 
law, but then decreases as the plunging trajectory nears 
the event horizon. Because we neglect higher multi- 
poles (which are extremely important in the fast-motion, 
strong-gravity region), this method likely underestimates 
the recoil. The total accumulated recoil at the cutoff time 
T using this method is shown in the lower curve of Figure 
(also for 7/ = 0.1). 

The shaded region between the two curves in Figure[2l3 
represents our uncertainty in the total recoil at the end of 
the plunge. This uncertainty is largest for retrograde or- 
bits around rapidly spinning holes, in which the distance 
the particle must "plunge" is greatest. For prograde in- 
spiral into rapidly spinning holes, much of the recoil is 
due to emission during the slow inspiral phase, for which 
our BH perturbation techniques are well-suited. Figure 
^ shows the relative contributions from the inspiral and 
plunge for such a scenario. 

Although the two calculations for the plunge recoil give 
rather different results, useful astrophysical information 
is contained in the approximate upper and lower bounds 
that they represent. The estimate V ~ 120km/s bisects 
the shaded region of Figure |2b and represents a typical 
recoil velocity for this mass ratio. Note also that the 
numbers in Figure El can be scaled to higher mass ratios 
by multiplying by f{q)/f{f] ~ 0.1). For q « 0.38 this 
implies that our results can be augmented by a factor 
w 2.3. 

5. DISCUSSION 

The punchline of this analysis is simple: quasi- 
Newtonian estimates have significantly overestimated the 
kick velocity from anisotropic GW emission during bi- 
nary coalescence. The recoil is strongest when the 
smaller member is deep in the strong-field of the large 
black hole. General relativistic effects, such as the gravi- 
tational redshift and spacetime curvature-scattering, act 
on the emitted GWs and reduce the recoil. 

Though reduced, the recoil remains large enough to 
have important astrophysical consequences. Recoils with 
V ~ 10-100 km/s are likely; kicks of a few hundred km/s 
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are not unexpected; and the largest possible recoils are 
probably < 500 km/s. These speeds are smaller than 
most galactic escape velocities, suggesting that BH merg- 
ers that follow galaxy mergers will remain within their 
host structures. However, these recoils are similar to the 
escape speeds of dwarf galaxies; and they may be suffi- 
cient to escape from mergers in high re dshift structures 
[z>5 - 10; cf. iBarkana fc Loebl l)200l|) . Fig. 8]. Binary 
BH ejection from globular clusters is quite likely, with 
significant implications for the formation of intermedi- 
ate mass black ho les (IMBH) via hierarchical mergers 
tMil ler fc Colbert! 120031. Our recoil estimates will also 
be useful in simulations of supermas sive and IMBH evo- 
lution in dark halos (Volonteri. Haar dt. fc Madau1l2003l: 
iMadau et al. "20041' 

Future papers will present the formalism used for this 
analysis, and will investigate the influence of orbital in- 
clination on the recoil. More work in perturbation theory 
also remains in addressing the recoil from the plunge and 
final ringdown of the merging black holes. 

Finally, Rcdmount & Rees ( 1989) have speculated that 
spin-orbit misalignment could lead to recoil directed out 
of the orbital plane. This recoil might accumulate sec- 
ularly rather than oscillate, and would be similar to 
the "electromagnetic rocket" in pulsars with off-centered 
magnetic dipole mom ents fHarrison & Tademaru 1973 
iLai. Chernoff. fc Co rdes 2001). We suspect that this ef- 
fect occurs but it is likely small compared to the recoil 
from the final plunge and merger. Firm estimates of the 
final kick velocity will rely on correctly modelling the 
final phase of BH coalescence. For comparable mass bi- 
naries, full numerical relativity will ultimately be needed 
to accurately compute the GW recoil. 
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